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Part  1  -  The  theory  of  stick-slip  is  developed,  based  on  the  concept  of  a 
static  and  dynamic  coefficient  of  friction.  The  dynamic  coefficient  is  as¬ 
sumed  to  be  independent  of  displacement  and  to  depend  wholly  on  normal  stress 
across  the  surface.  It  is  also  assumed  to  be  velocity-independent.  The  the¬ 
ory  predicts  that  the  stress  drop  during  stick-slip  is  independent  of  machine 
stiffness.  The  displacement  during  stick-slip  is  by  contrast  directly  pro¬ 
portional  to  the  machine  compliance.  Two  series  of  experiments  with  different 
fault  angles  would  serve  to  evaluate  the  dynamic  coefficient  of  friction. 

Part  II  -  Frictional  sliding  on  sawcuts  and  faults  in  laboratory  samples  of 
various  silicate  rocks  is  markedly  temperature-dependent.  At  pressures  from 
1  to  5  kb,  stick-slip  gave  way  to  stable  sliding  as  temperature  was  increased 
200®  to  500“  C.  The  particular  temperature  of  transition  to  stable  sliding 
varied  with  rock  type. 
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INTRODUCTION 

This  report  describes  our  t?chnical  findings  and 
accomplishments  during  the  first  six-month  period  of  the 
contract.  The  report  is  in  two  parts,  which  may  eventually 
be  published  separately  and  are  therefore  treated  as  if 
they  were  independent.  In  Part  I,  the  results  of  our 
analytical  studies  are  presented.  In  Part  II,  we  give 
the  results  of  our  measurements  of  friction  at  high 
temperature. 
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PART  I 

An  analysis  of  stick-slip  on  rock  surfaces  in 
the  laboratory 

Pierre-Yves  F.  Robin 


ABSTRACT 

The  theory  of  stick-slip  is  developed,  based  on  the 
concept  of  a  static  and  dynamic  coefficient  of  friction. 

The  dynamic  coefficient  is  assumed  to  be  independent  of 
displacement  and  to  depend  wholly  on  normal  stress  across 
the  surface.  It  is  also  assumed  to  be  velocity- independent . 
The  theory  predicts  that  the  stress  drop  during  stick-slip 
is  independent  of  machine  stiffness.  The  displacement 
during  stick-slip  is  by  contrast  directly  proportional 
to  the  machine  compliance.  Two  series  of  experiments  with 
different  fault  angles  would  serve  to  evaluate  the  dynamic 
coefficient  of  friction. 

INTRODUCTION 

Stick-slip  is  a  type  of  relative  motion  between  two 
surfaces  in  contact,  under  the  action  of  shear  stresses. 

It  is  characterized  by  jerky  movements,  separated  by 
intervals  during  which  no  significant  relative  displacement 
occurs.  This  effect  has  been  studied  in  metals  (1,  2,  3) 
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and,  more  recently,  in  rocks  (4,  5,  6).  Stick-slip  in  metals 
has  usually  been  interpreted  in  terms  of  two  coefficients  of 
friction,  one  static  ar d  one  dynamic.  However,  because  the 
experimental  systems  used  were  dynamically  complex,  the 
problem  of  deducing  the  value  of  the  dynamic  coefficient 
of  friction  from  the  observed  motion  has  not  yet  been  solved 
(7).  On  rock  surfaces,  stick-slip  has  also  been  interpreted 
as  the  result  of  a  dynamic  coefficient  of  friction  smaller 
than  a  static  one  (8) .  It  is  also  possible,  however,  that 
stick-slip  results  from  random  variations  of  the  average 
coefficient  of  friction  as  the  two  surfaces  move  by  each 
other. 

An  understanding  of  this  phenomenon  is  in.portant  to 
the  geophysicist,  because  stick-slip  is  a  possible  mechanism 
for  earthquakes  (6,  9,  10).  As  Eurridge  and  Knopoff  (10) 
forcefully  put  it,  "the  nature  of  the  friction  during  a 
shock  determines  the  configuration  of  the  system  when  it 
has  finally  come  to  rest.  It  is  this  final  state  that 
determines  the  conditions  surrounding  the  next  succeeding 
shock.  Hence,  if  the  demonstrations  of  the  laboratory  end 
numerical  models  are  borne  out  in  nature,  it  would  seem 
likely  that  the  nature  of  the  friction  on  a  fault  surface 
determines  the  statistical  properties  of  the  earthquake 
shocks  that  are  observed  . . . . " 

This  paper  analyses  the  dynamics  of  stick-slip  for 
the  experimental  system  which  is  most  likely  to  provide  data 
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pertinent  to  earthquake  studies.  The  analysis  shows  how 
characteristics  of  friction  may  be  obtained  from  experimental 
results.  Available  data  (11)  do  not  contradict  the  hypotheses 
proposed. 


REPRESENTATION  AND  CHARACTERISTICS  OF  STICK-SLIP  MOTION 

In  a  typical  experimental  study  of  friction  on  rocks, 
a  cylindrical  specimen  has  a  pre-existing  cut  at  an  angle 
with  its  axis.  Jacketed,  this  specimen  is  submitted  to  a 
constant  confining  pressure  p,  and  a  load  parallel  to  the 
axis  of  the  cylinder  is  then  applied.  The  displacement 
measured  is  the  relative  displacement,  parallel  to  the  axis, 
of  two  parts  of  the  machine  on  either  side  of  the  specimen, 
usually  sufficiently  far  from  the  actual  cut  to  show  no 
significant  displacement  during  the  stress  drop.  A  typical 
record  is  shown  on  Figure  1(a).  Such  an  experiment  is  often 
reported  as  Figure  1(b)  (e.g.  6),  obtained  from  Figure  1(a) 
by  removing  from  the  displacement  the  part  which  is  due  to 
the  elastic  compliance  of  the  machine.  Slopes  of  lines  like 
F1I2  are  the  same  as  the  slope  corresponding  to  the  elastic 
deformation  of  the  whole  specimen,  except,  perhaps,  very 
close  to  Ii.  It  is  concluded  that,  in  general,  no  significant 
sliding  occurs  along  the  cut  when  the  stress  increases  from 
Fi  to  l2*  Hoskins  et  al  (5),  with  an  experimental  system 
widely  different  from  the  one  described  here,  and  with  a 
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m  n  1 


Various  representations  of  a  stick-slip  experiment. 
V7e:;ter]y  Granite  (sample  diameter;  0.6?5  in;  ground 
surface,  at  angle  0  =  30°;  confining  pressure;  4.59  kb). 
After  Eyeries  (11)  . 
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normal  load  never  exceeding  70  bars,  showed  some  movement 
along  the  surface,  at  constant  value  of  stress,  before  the 
slip  occurs.  In  general,  however,  we  will  assume  that 
displacement  only  occurs  during  the  stress  drop. 

If  the  angle  of  cut  9  is  defined  as  shown  on  Figure  4, 
differential  stress  (Oi-p)  and  confining  pressure  p  can  be 
converted  into  normal  eind  shear  stresses  by  the  following 
equations 

=  p  +  (oi  -  p) .  Sin^e  (1) 

T  =  (.0 1  -  p)  .  Sin  6.  Cos  6  (2) 

Stick-slip  experiments  may  then  be  reported  as  on  Figure 
1(c),  where  the  displacement  plotted  is  the  relative 
displacement  a  parallel  to  the  surface,  or  its  axial 
equivalent  u, 

Au  =  Aa.  Cos  0  (3) 

Figure  1(c)  can  be  considered  as  transformed  from  Figure 
1(a)  by  subtracting  the  di'r.placement  due  to  the  compliance 
of  the  "machine",  where  "the  machine"  is  now  understood  to 
include  the  length  of  the  specimen.  Because,  indeed,  for 
the  surface  in  motion,  it  is,  to  this  point,  immaterial 
whether  the  compliance  is  provided  by  steel,  granite  or 
Teflon,  this  representation  is  preferable.  It  may  be 
pointed  out  also  that  the  affine  transformation  from  1(a) 
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to  1(c)  can  be  performed  geometrically,  without  having  to 
measure  separately  the  various  compliances  involved.  The 
displacement  of  the  driving  mechanism  during  the  time  of  a 
stress  drop  can  be  neglected,  and  the  slope  of  lines  like 
IiFi  (in  axial  load  vs.  axial  displacement  coordinates)  is 
therefore  equal  to  the  inverse  of  the  compliance  S  of  the 
"machine” . 

The  shear  stress  at  points  like  Ii  is  the  shear  stress 
necessary  to  overcome  friction.  After  a  number  of  stress 
drops  which  depends  on  the  confining  pressure,  some 
equilibrium  is  reached,  and  the  stresses  at  points  like 
I3,  1 1*,  etc.  remain  approximately  constant,  sometimes 
remarkaoly  so  (5,  6).  Values  of  the  stress  drops  are  often 
quite  regular;  they  are,  typically,  of  the  order  of  half  the 
maximum  stress. 

THEORY 

In  this  paragraph  and  the  next,  it  is  often  easier  to 
deal  with  differential  axial  load  F  and  axial  displacement  u 
than  with  the  shear  stress  t  and  the  displacement  parallel 
to  the  friction  surface  a.  As  seen  on  Figure  1(c),  the 
transformation  involves  only  a  change  of  scale  along  the 
axes,  and  the  conversion  factors  are  based  entirely  on  the 
geometry  of  the  system.  In  particular  u  (Eq.  3)  designates 
the  axial  component  of  a,  i.e.  of  the  relative  displacement 
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of  two  points  immediately  adjoining  the  cut  and  on  either 
side  of  it.  Also,  the  following  quantities,  having  the 
dimension  of  a  work,  are  equal: 

F.du  =  A.T.da 

where  A  is  the  surface  area  of  contact. 

Clearly,  stick-slip  requires  that  the  resistance  to 
shear,  immediately  after  initiation  of  motion,  wou].d  be 
smaller  than  before  motion.  If  this  were  not  the  case,  as 
an  infinitesimal  motion,  du,  causes  a  decrease  dF  =  S.du  in 
the  force  applied  by  the  machine,  further  motion  would  then 
be  resisted  by  friction,  until  the  driving  mechanism  raises 
F  again.  Sliding  on  the  surface  would  then  be  stable. 

Before  exploring  the  various  possibilities,  it  is  useful 
to  study  the  energy  transformations  during  stick-slip.  The 
argument  given  here  is  similar  to  the  one  given  by  Rabinowicz 
(7) ,  but  includes  the  seismic  energy  in  the  system. 

Let  Ej  be  the  initial  elastic  stored  in  the  machine, 
before  slip  occurs.  During  the  stress  drop,  neglecting  the 
very  small  amount  of  energy  given  by  the  driving  mechanism, 
this  energy  is  transformed  into  several  terms: 

El  =  Ej,  +  Ep  +  +  Eg  (4) 

where  is  the  heat,  or  friction  energy,  generated  on  the 
surface ; 
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Ep  is  the  sum  of  the  potential  elastic  energies  of  all 
parts  of  the  system; 

is  the  sum  of  the  kinetic  energies  of  all  parts  of 

the  system;  E  may  include  attenuation  energy,  i.e.  heat 
i\ 

generated,  by  vibrations,  elsewhere  than  on  the  friction 
surface; 

E^  is  the  work  done  against  the  confining  pressure. 

In  general,  when  relative  motion  across  the  friction 
surface  stops,  not  all  parts  of  the  system  stop.  The  terms 
Ep  and  of  equation  (4)  can  be  redistributed  into  an 
elastic  potential  term  for  the  final  position  E  ,  and  a 
seismic  term  Eg.  Eg,  like  Ej^,  may  include  or  transform  into 
heat,  generated  elsewhere  than  on  the  friction  surface.  The 
energy  equation  becomes 


El  =  E„  +  +  Eg  +  Ej, 


(5) 


Now,  the  origin  of  the  axial  displacement  u  can  be  taken  at 
the  position  before  motion,  and  u  may  be  defined  as  positive 
as  in  Figure  1.  The  heat  generated  is  then  given  by 


.  dy 


(6) 


Because  forces  are  taken  as  positive  when  corresponding  to 
compressions,  the  compliance  S  of  the  mv^chine  is  a  negative 
quantity.  It  can  be  assumed  constant  over  the  range  of  the 
stress  drop.  Calling  the  axial  force  on  the  surface  at 
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equilibrium,  before  initiation  of  motion,  the  decrease  in 
elastic  energy  is 

*  I  <’> 

o 

Substituting  (6)  and  (7)  into  (5)  gives: 


F)  dy 


(8) 


E  ,  the  energy  dissipated  seismically,  is  always  positive. 
Therefore  Equation  (8)  is  expressed  by  the  following 
inequality  for  the  oriented  surface  areas  of  Figure  2 ; 


Aj  +  Aj  >  0  (8' ) 

The  possibility,  mentioned  earlier,  that  stick-slip 
would  result  from  random  variations  of  the  average  friction 
coefficient  across  the  interface  can  be  represented  as  on 
Figure  3(a).  Several  features  of  Figure  3(a)  distinguish 
it  from  Figure  1(c)  and  make  this  theory  difficult  to 
reconcile  with  the  observed  phenomenon  without  further 
modification. 

After  a  stress  drop,  the  shear  stress  across  the 
surface  increases  without  displacement,  along  branches  like 
FiGi,  and  this  is  similar  to  branches  Fili  in  Figure  1(c). 

But  stable  sliding  along  branches  like  Fili  v/ould  then  occur, 
following  the  random  friction  curve.  The  stress  levels  at 


Energy  equation  for  stick-slip  when  friction  varies 
with  displacement. 


Irregular  stick-slip,  caused  by  variation  of  the  resistance  of 
friction  with  displacement. 
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point  I  or  F  would  also  be  expected  to  show  wider  variations 
than  the  ones  observed.  Finally,  a  reduction  of  the  compliance 
of  the  machine  would  be  expected  to  cause  stick-slip  to 
disappear,  ard  this  is  not  found  by  experiments  (6).  This 
model  can  account,  however,  in  a  rough  qualitative  way,  for 
the  observed  increase  in  the  amplitude  of  the  stress  drop 
with  cin  increase  in  confining  pressure.  Figure  3(b)  is 
similar  to  Figure  3(a),  drawn  for  the  same  hypothetical 
interface,  under  a  higher  normal  stress.  The  compliance 
of  the  machine,  however,  is  approximately  unchanged.  The 
stress  drops  can  be  seen  to  be  statistically  higher  than 
in  the  case  of  Figure  3(a).  In  spite  of  this,  and  because 
of  the  reasons  given  previously,  irregular  variation  with 
displacement  of  the  average  friction  as  an  explanation  of 
stick-slip  (Rabinowicz's  irregular  stick-slip,  7)  is 
tentatively  rejected  here.  The  other  possibility,  accepted 
for  metals,  is  to  consider  the  existence  of  a  dynamic 
coefficient  of  friction,  smaller  than  the  static  coefficient 
of  friction,  without  at  first  trying  to  explain  the  physical 
nature  of  this  difference.  A  priori,  many  factors  may 
influence  this  dynamic  coefficient  of  friction.  As 
suggested  by  its  very  appearance,  the  dynamic  coefficient 
may  depend  on  the  velocity  of  relative  displacement  during 
the  stress  drop  itself;  or  it  may  also  depend  on  the  amount 
of  relative  displacement  since  the  beginning  of  the  individual 
stress  drop.  The  dynamic  coefficient  may  also  depend  on  the 
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total  amount  of  displacement  since  the  beginning  of  the 
experiment.  It  may  be  a  function  of  the  normal  stress 
across  the  sliding  surface,  or  a  function  of  the  orientation 
of  the  stress  quadric  with  respect  to  the  surface. 

Before  developing  the  model,  for  reasons  of  simplicity, 
the  following  assumptions  will  be  made: 

(1)  The  dynamic  coefficient  of  friction  is  independent 
of  the  total  amount  of  displacement  since  the 
beginning  of  the  experiment.  This  may  be  justified 
by  the  fact  that,  after  a  few  initial  stress  drops, 
stick-slip  becomes  a  fairly  stable  phenomenon,  as 
mentioned  in  Section  II. 

(2)  The  dynamic  coefficient,  like  the  static  (e.g.  4, 

5,  12) ,  depends  only  on  the  normal  stress  across 
the  surface,  and  does  not  otherwise  depend  on  the 
orientation  of  the  stress  quadric.  The  law  of 
static  friction  for  many  rocks  can  be  expressed 
by  a  relation  of  the  form  (4,  5,  6,  12) 

^  *  "s  •  <9) 

where  S  and  y  are  constants.  In  the  same  way, 
s  s 

dynamic  friction  is  assumed  to  be  of  the  form 
^  ®d  *  °n 

is  constant,  i.e.  independent  of  the  amount 
of  displacement  during  the  individual  stress  drop. 


(3) 
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or  of  the  displacement  velocity.  It  is  not 
essential  that  M,  reaches  this  constant  value 

Cl. 

at  the  immediate  initiation  of  motion,  but  it 
is  assumed  tc  have  reached  it  during  most  of 
the  motion.  This  assumption  can  only  be 
justified,  refined  or  discarded  by  the  agreement 
or  disagreement  of  predictions  based  on  it  with 
experimental  results. 

The  developmei  t  of  the  model  is  now  a  relatively  simple 
matter. 

DEVELOPMENT  OF  THE  MODEL 

A  schematic  diagram  of  the  machine  of  the  model  is 
\  shown  on  Figure  4.  The  shear  stress  component  on  the  surface 
is  applied  through  a  "column"  of  springs  of  varying 
characteristics  (steel,  rock,  etc.).  The  surface  area  of 
friction.  A,  is  assumed  not  to  vary  significantly  during 
one  single  stress  drop.  Under  the  highest  confining  pressure, 
A  may  in  fact  vary  by  3  or  4%.  The  confining  pressure  p  is 
kept  constant  during  the  stress  drop. 

Immediately  before  initiation  of  motion  the  stress 
components  near  the  surface  are  related  by  Equations  (1)  and 
(2);  and  t  are  themselves  related  by  Equation  (9). 

Equations  (1),  (2)  and  (9)  combine  into 

'^*,s  P  ^  sin  6(Cos  0  -  M  Sin  6) 

s 


(11) 
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Consider  now  the  state  of  stress  during  motion,  in  the 
immediate  vicinity  of  the  surface.  The  inertial  forces  due 
to  lateral  accelerations  in  the  specimen  are  neglected.  The 
minimum  principal  stresses  are  still  0i  =  02  =  p.  The  angle  6 
being  unchanged,  Eqviations  (1)  and  (2)  are  still  the 
transformation  laws.  When  they  are  combined  with  Equation  (10), 
the  same  algebra  as  in  the  static  case  leads  to  a  similar 
equation: 

_  ^d  ^  ^d‘^ 

^i,d  ^  ^  Sin  6 (Cos  0  -  y^.Sin  6) 

To  0,  _  and  0  ,  correspond  axial  forces  F  ,  F, ,  by  the 

geometric  relation  F  =»  A  Sin  0.  oi.  The  problem  can  now  be 
expressed  in  the  following  manner:  a  spring  column,  along 
which  mass  and  compliance  are  distributed  non-uniformly ,  is 
compressed,  in  equilibrium,  by  a  force  F  .  The  force  is 
suddenly  reduced  to  a  lower  value  F^.  What  is  the  amplitude 
of  the  first  uniform  motion  of  the  end  of  the  spring? 

Figure  5  is  another  sketch  of  the  spring  column.  The 
coordinate  along  it  is  x.  The  column  is  fixed  rigidly  at 
X  =  0.  The  displacement  of  a  point,  u(x,  t) ,  is  a  function 
of  time.  The  specific  mass  per  unit  length  is  a  function  of 
position,  m(x) ,  and  is  always  positive.  The  specific  compliance 
is  also  a  function  of  position,  s(x).  Because  of  our  convention 
that  compressions  are  positive  and  that  displacements 
corresponding  to  stress  drops  are  also  positive,  s(x)  is 


Distribution  of  mass  and  compliance.  Case  of^a  rigid  mass  at 
the  end  of  a  massless  spring. 
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always  negative;  s (x)  can  be  zero  (rigid  section),  but 
cannot  be  infinite.  The  total  comftliance  of  the  machine 
is  given  by: 

J  s(x).dx  =  S,  negative  (13) 

o 


The  origin  of  displacements  u(x,  t)  may  be  talcen  as  the 
position  cf  each  section  immediately  before  slip.  It  can 
b*"'  shown  that  rhe  equation  of  motion  is: 


9^u  .1  9,1  9u. 

HTTx)  •  IST  ^il3E)  *  97^ 


0 


(14) 


The  initial  conditions  are  that 


at  t  =  0  : 


u(o,  t)  *  0 
9u  9^u 


9t 

1 


9t! 


=  0  for  any  x  ) 


9u 


(15) 


The  boundary  conditions,  when  t  >  0,  are 


u(o,  t)  =  0 


1  9u(l,t) 

s  (1)  ’  53c 


=  [F^  - 


F  ] ,  for  any  t  ) 

S  \ 


(16) 


Let  T  be  the  time  at  which  the  velocity  of  the  end  of 
the  spring  comes  to  0.  At  this  time,  the  coefficient  of 
friction  resumes  its  static  value  and  relative  displacement 
stops.  The  axial  component  of  the  relative  displacement 
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for  the  stress  drop  is  thus  u(l,  T) .  In  the  general  case. 
Equation  (14)  is  not  easy  to  solve  and  must  be  integrated 
numerically.  However,  if  the  mass  and  the  compliance  are 
distributed  as  in  Figure  6,  the  problem  becomes  one  of  an 
intertial  pendulum  with  a  massless  spring;  Equation  (14) 
may  be  solved  lirectly  (2,  7,  8).  Instead  it  can  also  be 
noted  that  for  such  a  distribution  of  mass  and  compliance 
no  seismic  energy  remains  in  the  system  after  the  end  of  the 
motion.  In  such  case.  Equation  (8)  gives  (see  also  Figure  7) : 

u(l,  t)  =  2(F^  -  Fg].S  (17) 

The  argument  leading  to  Equation  (8)  also  indicates  that 
this  is  the  maximum  value  of  u(l,  T)  for  a  machine  with  a 
given  total  compliance  S.  In  general 

u(l,  T)  =  K.S. [F^  -  Fg]  (17' ) 

where  0  <  K  <  2 

K  is  a  constant  coefficient  characteristic  of  the  machine, 
as  can  be  shown  easily.  If  uj  is  a  solution  of  Equation  (14) 
for  a  value  [F^  "  Boundary  Condition  (16)  ,  U2  =  a  ui 

is  a  solution  for  a  boundary  value  a[F^  “  particular, 

if  Ui  (1,  T)  ,  equal  to  Ki.S.lF^  -  ,  is  the  displacement  in 

the  first  case,  the  displacement . in  the  second  case  is; 

Uzdf  T)  =  a.uid,  T)  =  a.Ki.lF^  -  F^]  =  Kj.S  a(F^  -  F^]. 
Therefore  Kj  =  Kj  =  K  . 


Stick-slip  for  a  machine  having  a  distribution  of  mass  and 
compliance  as  in  Figure  6. 
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Increase  of  stress  drop  with  confining  pressure  data  from 
Byerlee  (11) .  Corrected  for  the  decrease  in  contact  area. 
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The  load  drop  is  what  is  recorded  directly  in  an 
experiment.  In  absolute  value,  it  is  equal  to: 


AL  =  KlF,  -  F^] 
a  s 


Using  Equations  (11)  and  (12) ,  it  is: 

S,  S 

AT.  K  A  f  ^  ® 

^Cos  e  -  y,.Sine  Cos  6  -  y  .Sin  0 

Q  S 

^  'Cos  6  -  y,.Sin  6  ”  Cos  6  -  y  .sin  6^ ^ 
a  s 


(18) 


The  load  drop  is  thus  found  to  be  independent  of  the 
compliance  S  of  the  machine.  This  was  found  experimenta) ly 
by  Byerlee  and  Brace  (6) .  Curves  obtained  by  Byerlee  (li) 
for  friction  of  Westerly  granite  (with  6  =  30°)  show  that 
AL  does  increase  linearly  with  confining  pressure  as 
predicted  by  Equation  (18). 

The  quantities  and  are  readily  obtained 
experimentally,  and  have  been  obtained  in  particular  for 
Westerly  granite  (12)  .  The  only  unlcnowns  in  Equation  (18) 
are  therefore  K,  and  y^.  The  slope  r  of  the  line  in 
Figure  8  gives  the  following  equation; 

Pd 


r  =  K  A( 


Cos  6  -  y^.Sin  6 


^  \ 
Cos  6  -  y  .Sin  0' 


(19) 


or,  numerically: 


1.28(cm2)  =  7.92  K  (p  ~p - 0.53)  (cm^)  (19') 
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Another  series  of  experiments,  with  a  different  angle 
0',  would  give  another  equation  similar  to  (19).  The  two 
equations  would  then  be  sufficient  to  determine  K  and 
The  intersection  of  the  two  line3  (as  in  Figure  8)  with  the 
ordinate  axis,  p  =  0,  should  both  give  the  same  value  of  S^. 
There  is,  however,  a  large  degree  of  uncertainty  for  these 
intersections. 
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PART  II 

Experimental  studies  of  high  temperature  friction 
W.F.  Brace  and  D.K.  Riley 


ABSTRACT 

Frictional  sliding  on  sawcuts  and  faults  in 
laboratory  samples  of  various  silicate  rocks  is  markedly 
temperature-dependent.  At  pressures  from  1  to  5  kb, 
stick-slip  gave  way  to  stable  sliding  as  temperature  was 
increased  200®  to  500®C.  The  particular  temperature  of 
transition  to  stable  sliding  varied  with  rock  type. 


Several  field  and  laboratory  observations  suggest 
that  earthquakes  result  from  a  large-scale  form  of  stick- 
slip  (5) .  For  one  thing,  unstable  (stick-slip  and  earth¬ 
quake-producing)  and  stable  (fault  creep  and  stable  sliding) 
motion  have  been  found  both  in  the  field  and  in  laboratory 
experiments  on  rocks  under  high  pressure.  For  another,  the 
same  mineralogic  controls  on  stability  have  been  noted  in 
field  euid  laboratory.  For  example,  minute  amounts  of 
serpentine  in  a  dunite  produced  stable  sliding  (6) ;  fault 
creep  in  California  seems  restricted  to  areas  whero  the  San 
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field  euid  laboratory.  For  example,  minute  amounts  of 
serpentine  in  a  dunite  produced  stable  sliding  (6) ;  fault 
creep  in  California  seems  restricted  to  areas  where  the  San 
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Andreas  fault  system  cuts  serpentine-bearing  rocks  of  the 
Franciscan  series  (7) . 

How  can  the  disappearance  of  earthquakes  at  shallow 
depths  be  explained  on  the  basis  of  laboratory  studies? 

Three  possibilities  are  apparent:  a  mineralogic  change 
with  depth,  existence  at  depth  of  certain  pore  pressure 
conditions  known  to  stabilize  sliding  in  rocks  (8) ,  and 
temperature  increase.  The  last  is  the  least  understood. 
Stable  sliding  at  high  temperature  is  suggested  by  a  few 
observations  of  stable  faulting  at  high  temperature  (9)  and 
by  somewhat  ambiguous  results  with  powders  deformed  between 
rotating  emvils  (10) . 

There  are  relatively  few  laboratory  studies  of  rock 
fracture  at  high  temperature  and  pressure  (11)  and 
practically  none  of  frictional  sliding.  The  biggest 
experimental  difficulty,  particularly  for  sliding,  is 
jacket  design.  The  jacket,  required  in  a  triaxial  experiment 
to  exclude  the  gas  pressure  medium  from  the  rock  sample,  is 
typically  metal  foil.  The  foil  ruptures  easily  at  any 
strain  discontinuity  such  is  a  fault.  Our  procedure  was 
to  retain  the  thin  foil,  but  to  add  a  sleeve  of  graphite 
between  rock  sample  and  foil  (Figure  1) .  The  sharp  offset 
at  the  fault  is  smaared  out  in  the  soft  gr  phite,  and 
appreciable  motion  on  the  fault  is  tolerated  before  the 
foil  ruptures.  This  simple  modification  in  jacket  design 
enabled  us  to  use  otherwise  standard  experimental  procedure 
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for  high  temperature  deformation  study,  to  obtain  the  first 
detailed  picture  of  the  effects  of  temperature  on  friction 
of  rocks . 

The  rock  sample  was  a  precisely  ground  cylinder  16  mm 
in  diameter  by  35  mm  long.  Sawcut,  if  present,  was  located 
midway  between  the  ends  and  made  an  angle  of  30®  to  the 
cylinder  axis.  Graphite  sleeve  was  1.3  mm  thick,  the 
annealed  seamless  copper  foil,  0.32  mm  thick. 

An  extensive  series  of  experiments  was  conducted  at 
room  temperature  to  determine  any  possible  stabilizing 
effect  of  the  graphite-copper  jacket.  The  results  (Figure  2) 
revealed  that  a  stabilizing  effect  on  sliding  existed  only 
below  2  kb  pressure.  This  took  the  form  of  lowering  the 
amplitude  of  stick-slip  to  nearly  zero.  The  shearing  stress 
to  cause  frictional  sliding  increased  about  10  percent 
compared  with  an  experiment  at  the  same  pressure,  using  a 
3  mm  thick  polyurethane  jacket.  Because  of  these  effects, 
most  experiments  here  were  conducted  at  or  above  2  kb  pressure 
where  room  temperature  experiments  using  copper-graphite  were 
nearly  identical  with  those  using  polyurethane.  In  any 
event,  stabilizing  effects  present  at  room  temperature 
probably  would  not  be  important  at  high  temperature  because 
of  increased  ductility  of  the  copper. 

Our  apparatus  resembled  in  a  general  way  that  described 
on  page  46  of  Reference  9.  It  was  internally  heated,  with  low 
friction  0-ring  seals.  Stiffness  of  the  loading  system  was 
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about  10*  kg/cm.  Pressure  was  known  to  1  percent, 
temperature  to  about  10®.  Strain  rate  was  calculated  from 
the  rate  of  advance  of  the  screw-driven  piston. 

Frictional  sliding  was  studied  in  triaxial  experiments 
in  which  the  cylindrical  sample  contained  a  fault  or  sawcut 
(12) .  The  fault  was  formed  by  loading  an  initially  intact 
sample  to  failure;  the  sawcut  was  m!ade  in  the  sample  at  an 
angle  (30®)  close  to  that  of  typical  faults  (26-32®). 

Of  our  two  types  of  experiment,  with  sawcut  and  with 
fault,  presumably  the  latter  more  nearly  resembles  actv.al 
faults.  Sawcuts  are  flat  and  have  a  finely  ground  surface; 
faults  have  abundant  gouge  and  the  surface  irregularity  one 
normally  associates  with  actual  faults.  Unfortunately  a 
'fault'  experiment  is  more  difficult  and  the  results  often 
more  ^unbiguous  than  a  'sawcut'  experiment.  For  example, 
each  laboratory  fault  differs  in  detail;  sawcuts  are  nearly 
identical;  as  a  result,  data  from  faulted  samples  shows 
greater  scatter  than  that  from  sawcuts.  For  exploratory 
work,  results  from  sawcuts  are  probably  valid;  Byerlee  (13) 
found  for  granite  only  minor  differences  in  friction  between 
sawcuts  and  faults  once  some  motion  had  occurred.  Most  of 
the  results  reported  here  are  for  sawcuts. 

We  studied  frictional  sliding  in  Westerly  granite  and 
San  Marcos  gabbro  (14)  at  pressures  to  5  kb,  temperatures  to 
525®C,  strain  rates  from  10““  to  10“*sec”‘.  The  samples 
were  vented  to  the  atmosphere  through  a  hollow  piston,  so 
that  presumably  pore  pressure  was  nearly  zero. 
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The  results  for  both  rocks  are  shown  in  Table  1. 

In  the  first  column,  W  refers  to  the  granite  and  SM  refers 
to  gabbro.  The  displacement  xn  the  fourth  column  is 
displacement  at  pressure  and  temperature  on  the  fault. 

In  all  runs  the  sample  was  heated  at  temperature  for  one 
hour,  and  the  strain  rate  was  10“ ®  per  second  unless 
otherwise  noted. 

Results  are  shown  in  Figure  2  for  sawcuts  in  granite. 
Apparently  high  temperature  had  a  strong  stabilizing  effect 
on  stick-slip;  large  cimplitude  stick-slip  at  low  temperature 
(the  22®  curve  in  Figure  3)  gave  way  to  stable  sliding  as 
temperature  is  increased  (the  306®  curve  in  Figure  3) , 
Results  at  higher  pressure  were  similar.  No  change  in 
character  of  the  sliding  was  evident  over  the  3  mm  or  so 
of  sliding  motion,  which  was  the  limit  imposed  by  the 
apparatus.  Neither  strain  rate  nor  heating  procedure 
appeared  to  affect  behavior  cuch  as  shown  in  Figure  3. 
Samples  were  run  at  10”“  to  10”^ sec” ^  strain  rate,  were 
heated  at  pressure  for  1  to  25  hours,  with  and  without 
vacuum  (10“^  Torr) ,  and  were  heated  and  then  run  at  room 
temperature . 

Results  for  faults  in  gabbro  are  shown  in  Figure  4; 
these  are  also  typical  for  granite.  The  faults  were  formed 
in  the  samples  at  0.5  to  1  kb  pressure  and  room  temperature. 
Pressure  and  temperature  were  then  raised  to  the  conditions 
of  the  friction  experiment.  The  marked  effect  of  high 
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TABLE  1 


High  Temperature  Friction 


Expt . 

P 

kb 

T 

“C 

Displ . 

mm 

kb 

Motion 

Remar]  s 

W-l-S 

2.08 

470 

2.75 

2.65 

STA 

e  =  10  “/sec 

W-2-S 

2.09 

303 

2.40 

4.6 

STA 

e  =  10”‘*/sec 

W-4-S 

4.00 

293 

1.50 

8.5 

STK 

Vac.htg.l2  h  @  140“C 

W-5-S 

1.00 

303 

2.1 

3.2 

STA 

e  ranged  from  10"“ 

W-6-S 

4.03 

438 

1.4 

4.8 

STA 

to  10"®/sec 

W-7-S 

1.99 

305 

1.4 

5.3 

STA 

25  hr.htg.at  2  kb, 

W-7'-S 

3.05 

305 

0.  8 

7.9 

STA 

300“C 

W-9-S 

2.00 

191 

1.2 

5.5 

STA 

W-9'-S 

4.00 

197 

1.0 

9.1 

STK 

W-IO-F 

4.04 

439 

1.5 

11.3 

STA 

W-ll-F 

4.00 

260 

3.7 

11.7 

STK 

e  =  10~‘*/sec 

SM-12-F 

4.00 

398 

1.5 

11.0 

STA 

e  =  lO'^/secM  hrs. 

SM-13-F 

5.00 

197 

4.3 

11.5 

STK 

htg.  @  400®C 

SM-14-F 

4.00 

140 

2.70 

11.0 

STK 

e  =  lO'^/sec 

W-15-F 

1.92 

527 

2.30 

6.5 

STA 

W-16-5 

1.91 

140 

1.0 

3.  0 

STK 

W-16'-5 

1.91 

203 

0.8 

3.7 

STK 

SM-17-F 

4.04 

512 

2.8 

9.3 

STA 

SM-18-F 

4.13 

100 

2.0 

11.5 

STK 

VJ-19-F 

5.05 

403 

3.2 

13.9 

STA 

W-20-F 

6.11 

405 

1.6 

14.9 

STA 

W-21-F 

4.95 

350 

1.5 

12.3 

STK 

W-22-F 

3.00 

390 

1.5 

7.  8 

STA 

W-25-F 

3.52 

206 

2.2 

11.5 

STK 

W-26-F 

5.02 

512 

2.0 

13.1 

STA 

SM-27-F 

3.50 

355 

2.6 

8.3 

STA 

SM-28-F 

3.00 

103 

2.0 

9.0 

STK 

SM-29-F 

4.50 

306 

3.1 

9.6 

STA 

Imm 
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temperature  is  again  evident,  although  the  transition 
between  stick-slip  and  stable  sliding  appears  less  sharp 
than  for  sawcuts.  In  other  words,  there  appears  a 
significant  range  of  temperature  over  which  stick-slip  was 
preceded  by  some  stable  sliding.  It  is  not  certain  that 
even  at  the  highest  temperature  stick-slip  would  not  have 
occurred  had  there  been  additional  displacement. 

Results  for  both  sawcuts  and  faults  are  shown  in 
Figure  5,  in  which  stick-slip  (open  figure)  or  stable 
sliding  (closed  figure)  is  indicated.  Where  appreciable 
stable  sliding  preceded  the  stick-slip,  this  is  designated 
by  a  half-closed  figure. 

Several  features  are  evident  in  Figure  5.  First, 
sliding  on  granite  sawcuts  has  a  well-defined  field  of 
stability;  thus,  the  sliding  was  stable  at  high  temperature 
and  low  pressure,  and  unstable  at  high  pressure  and  low 
temperature.  Second,  the  field  boundary  for  the  sawcuts 
is  very  sharp;  within  about  lOC’C  large  amplitude  stick- 
slip  gave  way  to  stable  sliding.  Third,  the  results  for 
the  faults  in  granite,  although  very  limited  in  number, 
are  at  least  consistent  with  results  for  the  sawcuts;  the 
transition  from  unstable  to  stable  may  be  more  gradual  for 
faults  than  sawcuts  at  the  pressures  of  these  tests. 
Finally,  the  field  boundary  may  be  different  for  gabbro 
and  granite  with  stick-slip  disappearing  at  lower 
temperature  for  gabbro. 


F(G.  S 
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At  present  no  physical  explanation  can  be  offered 
for  this  pronounced  effect  of  temperature  on  s tick-slip. 

Our  understanding  of  the  stick-slip  process  is  still 
rather  incomplete.  Byerlee's  studies  (15)  suggest  t  lat 
brittle  fracture  plays  an  important  role  in  the  frictional 
behavior  of  rocks;  perhaps  an  explanation  would  be  a  ^parent 
if  more  were  known  about  effect  of  temperature  on  brittle 
fracture  of  rock-forming  minerals.  A  significant 
observation  from  the  present  work  is  that  frictional 
strength  is  lowered  by  temperature  by  abaut  the  same  amount 
as  fracture  strength,  relative  to  room  temperature  values. 
This  suggests  that  the  behavior  on  a  small  scale  is  the 
same  in  both  cases.  The  details  of  this  behavior  are  still 
obscure. 

Before  we  apply  present  results  to  real  faults,  we 
need  to  consider  differences  which  still  exist  betwe^in  the 
laboratory  experiment  and  the  field,  other  than  the  (Obvious 
one  of  scale.  Our  experiments  will  need  to  be  repeated 
with  pore  water  pressure,  for  presumably  natural  rocks  are 
wet.  Probably  the  effective  stress  law  will  be  followed 
as  it  is  at  room  temperature  (16).  Additional,  chem  .cal 
effects  may  be  present,  though,  to  judge  from  the  water¬ 
weakening  observed  in  certain  silicate  minerals  (17). 
Presumably  these  will  have  a  further  stabilii'ing  inf  uence 
emd  shift  the  field  boundary  of  Figure  3  to  tne  left,  to 
lower  temperatures.  Slower  strain  rate  than  used  heje  need 
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also  be  considered;  at  high  temperature  some  vvelding  or 
sintering  may  be  expected  at  very  slow  strain  rates,  and 
this  could  lead  to  stick-slip.  Further  study  is  needed 
here,  as  well.  Finally,  the  effects  of  displacement  will 
have  to  be  excunj.ned  more  fully.  To  judge  from  our  ooservation; 
with  faults  in  gabbro  and  granite,  the  nature  of  the  sliding 
motion  changes  somewhat  with  displacement  (the  middle  curve 
of  Figure  4).  Some  way  of  obtaining  much  larger  disj^lacements 
in  our  ladDoratory  seimples  is  needed. 
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